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ABSTRACT 



Despite the detection of a large number of gravitational microlensing events, 
the nature of Galactic dark matter remains very uncertain. This uncertainty 
is due to two major reasons: the lens parameter degeneracy in the measured 
Einstein timescale and the blending problem in dense field photometry Recently 
consideration has been given to routine astrometric followup observations of 
lensing events using the Space Interferometry Mission (SIM) as a means of 
breaking the lens parameter degeneracy in microlensing events. In this paper, 
we show that in addition to breaking the lens parameter degeneracy, SIM 
observations can also be used to correct for nearly all types of blending. 
Therefore, by resolving both the problems of the lens parameter degeneracy 
and blending, SIM observations of gravitational lensing events will significantly 
better constrain the nature of Galactic dark matter. 
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1. Introduction 

Surveys to detect Galactic dark matter by monitoring light variations of source stars 
caused by gravitational microlensing are being carried out by several groups and ~ 300 
events have been detected (MACHO: Alcock et al. 1997a, 1997b; EROS: Ansari et al. 1996; 
OGLE: Udalski et al. 1997). The light curve of a lensing event is represented by 
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where the lensing parameters f3, t , and t& represent the lens-source impact parameter, 
the time of maximum amplification, and the Einstein ring radius crossing time (Einstein 
timescale), respectively. Once the light curve of an event is observed, these lensing 
parameters are obtained by fitting the theoretical light curves of equation (1.1) to the 
observations. One can obtain information about individual lenses because the Einstein 
timescale is related to the physical lens parameters by 

r E / AGM DoiDis\ 1/2 
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where r E is the Einstein ring radius, v is the lens-source transverse speed, M is the lens 
mass, and D i, Di s , and D os are the separations between the observer, lens, and source star. 
However, since the Einstein timescale depends on a combination of the lens parameters, the 
values of the lens parameters determined from it suffer from large uncertainties. 

The lens parameters suffer from additional uncertainties due to blending. The 
probability for a single source star to be in the state of gravitational amplification, i.e., the 
optical depth r, is very low; from (O)l0~ 7 to (O)10~ 6 depending on the observed field. 
To increase the event rate, therefore, experiments are being conducted toward very dense 
star fields such as the Galactic bulge and Magellanic Clouds. However, photometry of such 
dense star fields is affected by the flux from stars not participating in the gravitational 
lensing, the "blending problem." Depending on the origin of the blended light, blending is 
classified into regular, amplification-bias, binary-source, and bright-lens blending (see § 2). 
When an event is affected by blending, the observed light curve is represented by 

A ohs = A ahs (l -/) + /; / = (1-3) 

-to + ts 

where Fq is the unblended flux of the source star before (or after) gravitational amplification, 
B is the amount of blended flux, and thus / represents the fraction of the blended light 
to the total observed flux. Therefore, when fitting a blended lensing event light curve, 
an additional lensing parameter / must be included. As a result, the determined lens 
parameters become even more uncertain. 

Efforts to improve the results of gravitational microlensing experiments are therefore 
focused on breaking the lens parameter degeneracy and correcting for the blending effect. 
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Various methods have been proposed to do this. The methods for removing the degeneracy 
in the lens parameters are summarized by Jeong, Han, & Park (1999). Here we discuss 
the methods proposed to correct for various types of blending (see § 2). However, these 
methods are either applicable to only in a few rare cases or impractical due to the large 
uncertainty in the determined value of /. Therefore, to gain improved results from the next 
generation of lensing experiments, it is essential to devise practical methods that can resolve 
the problems of the lens parameter degeneracy and blending in general microlensing events. 

Recently, routine astrometric followup observations of lensing events with high precision 
instruments such as the Space Interferometry Mission (hereafter SIM, Allen, Shao, & 
Peterson 1998) are being discussed as a method to break the lens parameter degeneracy 
of general microlensing events. When a source star is gravitationally lensed, it is split 
into two images. The image separation is too small for direct observation. However, the 
displacements of the light centroid of the two images can be measured with SIM. The 
trajectory of the centroid shifts is an ellipse (astrometric ellipse) whose shape depends on 
the lens-source impact parameter (see § 4.2). The usefulness of astrometric measurements 
of lensing events is that one can determine the angular Einstein ring radius, 9 E = r E /D Q i, 
from the measured centroid shifts because the size of the astrometric ellipse is directly 
proportional to 6> E . While the Einstein timescale depends on three lens parameters (M, 
D os , and v), the angular Einstein ring radius depends only on two parameters (M and D i). 
Therefore, by measuring 6 E , the uncertainties in the determined lens parameters can be 
significantly reduced (Walker 1995; Paczyhski 1998; Boden, Shao, & Van Buren 1998; Han 
k Chang 1998). 

In this paper, we show that in addition to breaking the lens parameter degeneracy, 
observations with SIM can be used to correct for the blending effect in general microlensing 
events. Due to the high resolution, <^ 0."5, from the diffraction-limited space observations, 
a significant fraction of events blended by field stars, i.e., regular and amplification-bias 
blended events, are automatically resolved by a single mirror of SIM. For events affected 
by very close blends with separations A9 <; 10 mas, including a large fraction of binary 
companions and all bright lenses, one can also identify the lensed source star by detecting 
the distortions of the astrometric centroid trajectory caused by the blended star. For events 
affected by medium range separation blends, 10 mas ^ A6 <^ 0."5, the blending effect can be 
corrected for by either detecting the shift of the broad envelope centroid between the point 
spread functions (PSFs) of lensed and blended stars caused by gravitational amplification 
or by examining the narrow fringe patterns observed at multi-wavelengths. Therefore, by 
resolving both the problems of lens parameter degeneracy and blending, SIM observations 
of gravitational lensing events will significantly better constrain the nature of Galactic dark 
matter. 
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2. Types of Blending 

Depending on the origin of the blended light, blending is classified into several types. 
First, "regular blending" occurs when a star brighter than the detection limit, which is 
set by crowding, is lensed and the measured flux is affected by the residual flux from 
other blended stars fainter than the detection limit. Due to regular blending, the apparent 
amplification of the event is lower than its intrinsic value. As a result, the apparent Einstein 
timescale is shorter than its true value, resulting in systematic underestimation of lens 
masses (Di Stefano & Esin 1995; Wozniak & Paczyhski 1997; Han 1998b). Since the optical 
depth is directly proportional to the timescale, the value of r determined without properly 
correcting for regular blending is also underestimated. 

The observed light curve is also affected by blending when one of several unresolved 
faint stars below the detection limit is lensed and its flux is associated with the flux from 
other field stars in the effective seeing disk of a monitored bright source star. This is 
known as "amplification-bias blending" (Bouquet 1993; Nemiroff 1994). The effects of 
amplification-bias blending for the determination of individual lens parameters are similar 
to those of regular blending, but the effects are more severe due to the much larger amount 
of blended light. However, the most important effect of amplification bias blending is that 
when not accounted for, the optical depth might be significantly overestimated (Alard 
1997). This is because r is determined based only on the number of monitored stars brighter 
than the detection limit, while events are actually detected among a larger number of stars, 
including fainter stars. Han (1997) estimated that the increase in r caused by the miscount 
of effectively monitored stars is large enough to compensate the decrease in r due to the 
decrease in timescales and make the observed optical depth overestimated by a factor ~ 1.7. 

In addition to background stars the lens itself can cause blending known as "bright 
lens blending" (Kamionkowski 1995; Buchalter & Kamionkowski 1997; Nemiroff 1997; 
Han 1998a). Bright lens blending causes the measured timescale and optical depth to 
be underestimated in the same manner as regular blending. In addition to the effects of 
regular blending, lens blending causes the determined optical depth to depend on the lens 
location. This is because detecting events caused by bright lenses close to the observer 
is comparatively more difficult than detecting events produced by lenses near the source. 
Additionally, since the brighter the lens is, the more massive it tends to be, and thus the 
more likely it is to be affected by lens blending. As a result, the decrease in optical depth 
for massive lenses is relatively bigger than the decrease for low-mass lenses, making the 
measured optical depth dependent on the lens mass function. 

Finally, the last type of blending, known as "binary star blending" occurs when the 
source is a binary. In general, the typical separation between the component stars in a 
binary system is larger than the typical size of the Einstein ring radius. In this case only 
one star is significantly amplified, and the flux from the other star simply contributes as 
blended light (Dominik 1998). For some binaries with small component separations, on the 
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other hand, both component stars can be amplified. However, even in these cases it is not 
easy to detect the binarity of the source because the light curve mimics that of a single 
source event with a longer timescale and larger impact parameter than the true values 
(Griest & Hu 1992; Han & Jeong 1998). 



3. Limitations of Current Blending-Correction Methods 

There have been various methods proposed to correct for the blending problem. 
However, these methods are either applicable to only a few specific types of blending or 
impractical due to the large uncertainty in the determined value of /. In this section, we 
list these proposed methods and discuss their performance in correcting blending effects. 

Because a blended event light curve is not exactly the same as that of an unblended 
event, ideally the effects are detectable using photometry. However, due to the limits on 
photometric precision of the current lensing experiments, blending has been photometrically 
detected only in a small number of events. Nevertheless, photometry is still being discussed 
because of the rapidly increasing photometric precision made possible through the operation 
of early warning systems (MACHO: Alcock et al. 1996; OGLE: Udalski et al. 1994) and 
subsequent effective followup observations (PLANET: Albrow et al. 1998; GMAN: Alcock 
et al. 1997c). However, we still find that even with this higher precision photometry, the 
uncertainties in the derived lensing parameters are very large, making it difficult to properly 
correct for blending effects. To demonstrate this, we simulate model events that are affected 
by various amounts of blended light and estimate the uncertainties of the recovered lensing 
parameters. The events are assumed to have (5 = 0.3, and £e,o = 15 days, which are the 
most common values for the Galactic bulge events, and are affected by blending with various 
blended light fractions of / = 0.3, 0.5, 0.7, and 0.9. The events are assumed to be observed 
5 times/day during — 0.5£e < ^obs < 3£e with a high photometric precision of p = 1%. The 
light curves are then fit with a theoretical blended light curve as given by equation (1.3). 
We estimate the uncertainties of the recovered lensing parameters by computing the values 
of x 2 by 



and theoretical model light curves, respectively. In Figure 1, we present the resulting 
values of x 2 as functions of £eAe,o and /. The degrees of freedom are determined by 
dof = Adat — A par — 1, where A par = 4 is the number of lensing parameters for the 
blended light curve fit. The contours are drawn at the levels of A% 2 = 1.0, 4.0, and 9.0 
(i.e., 1 a, 2 a, and 3 a levels) from inside to outside. From the figure, one finds that the 
uncertainties in both the derived values of £e and / are large. One also finds that the 
uncertainty ranges of the blended light fraction measured at all levels include / = (i.e., 




(3.1) 



where Ad a t is the number of data points and Aq and At represent the simulated 
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no blended light), implying that it will be difficult to detect blending effect even with high 
precision photometry. In addition, since the uncertainties increase rapidly with increasing 
values of /, correcting for blending effects using this method will be especially difficult for 
amplification-bias events which are severely affected by blending. 

Another method of correcting for blending effects is to detect color shifts during an 
event with high precision, multi-band photometry (Buchalter, Kamionkowski, & Rich 1996). 
This method might be applicable if the lensed source star has a very different color from 
that of the integrated light of the blended stars. However, the expected amount of color 
shift is generally very small because most Galactic bulge stars have similar colors (Goldberg 
1998). In addition, even when the color shift is detected, the determined lens parameters 
suffer from the same degeneracy as in the case of single-band photometry, making it difficult 
to properly correct for the blending effect (Wozniak & Paczyhski 1997). 

A more general method of correcting for blending effects is by detecting the shifts of 
a source star's image centroid, 5x, during an event (Alard, Mao, & Guibert 1995; Alard 
1996). If one of many stars in a blended seeing disk is gravitationally amplified, the position 
of the center of light will shift toward the lensed star by an amount 



where (x) is the position of the center of light before gravitational amplification and Xq is the 
location of the lensed star. Goldberg & Wozniak (1997) actually applied this method to the 
OGLE data base and found that 7 out of 15 tested events showed significant centroid shifts 
of 5x ^ 0".2, demonstrating the usefulness of this method. However, not all blended events 
produce large centroid shifts. If the amplification of an event is very low, i.e., ~ 1, the 
expected centroid shift is small since r] ~ 0. In addition, if the lensed star is the brightest 
one in the blended seeing disk and its flux dominates that from the other blended stars, the 
expected amount of centroid shift is very small even for a high amplification event because 
the position of the center of light before gravitational amplification will be very close to 
that of the lensed star, i.e., \{x) — xq\ ~ 0. For large centroid shifts, therefore, the lensed 
star should be one of the faint stars in the seeing disk so that it has a negligible effect on 
the position of the center of light (Han, Jeong, & Kim 1998). For amplification-biased 
events, source stars are generally very faint. Because they are faint, the fact that they 
are detectable implies that the source stars are highly amplified. Since the conditions for 
detecting amplification-biased events agree well with those for large centroid shifts, the 
centroid shift measurement is an efficient method for detecting amplification-bias blending. 
However, this method is not efficient for detecting other types of blending. First, because 
of a relatively small amount of blended light, regular blended events do not need to be 
highly amplified to be detected. Although they can be highly amplified, the dominance of 
the source flux over that from other faint blended stars will result in small centroid shifts. 
This method is also inefficient for an event affected by the light from very closely located 
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(3.2) 
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blends such as binary companions and bright lenses due to the small amount of expected 
centroid shift. 



In the previous section, we showed that correcting for blending effects for general 
microlensing events is difficult by using the previously proposed methods. In this section we 
show that with diffraction-limited, single-mirror images and the high astrometric precision 
of SIM, one can correct for the effects of nearly all types of blending. 

The optics of SIM can be modeled by a diffraction-limited telescope with two circular 
mirrors. Then the expected PSF is given by the classical double circular-aperture diffraction 
pattern formula of 



where Iq is the intensity of the central peak, d is the separation between the two mirrors 
with an aperture diameter D, A is the observed wavelength, 9 is the angle measured from 
the central fringe, and J\ is the lst-order Bessel function (Hecht 1998). According to the 
specification of SIM, d = 10 m and D = 0.3 m. The upper left two panels of Figure 2 show 
the expected PSF when a single (unblended) source star is observed by SIM. When observed 
at A ~ 5500 A, the PSF has a "broad envelope" of a width 1.22A/D ~ 0."5 modulated by 
"narrow fringes" of a width X/d ~ 10 mas. 

When the source star is blended, however, the PSF takes a different form from 
that in equation (4.1). The blended source star PSF can be classified into three regimes 
depending on the separation A0 between the lensed source and the blended star: regime 
1 (A0 £ 1.22A/D), regime 2 (A0 < X/d), and regime 3 (X/d < A6 < 1.22X/D). In the 
following subsections, we investigate how the blends with various separations affect the PSF 
and their effects can be corrected for individual cases. 



In this regime the broad envelopes of the two sources do not overlap (see the lower left 
panels of Figure 2), and thus the blended source is easily distinguished. In Baade's Window 
the density of stars with V < 23, which are the major sources of blending, is less than 
1 arcsec~ 2 according to the luminosity function determined from the Hubble Space Telescope 
observations by Holtzman et al. (1998). The stellar density in LMC field is even smaller. 
This implies that a significant fraction of sources blended by field stars, i.e., regular and 
amplification-bias blended events, are automatically resolved by a single mirror of SIM. 



4. Correction of Blending Effect by SIM 




(4.1) 



4.1. Regime 1 (A0 £ 1.22A/D ~ 0."5) 
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4.2. Regime 2 (A6> < X/d ~ 10 mas) 

In this case the narrow fringes from the two images overlap (the upper right panels 
of Figure 2). The position of the centroid of the central fringe, which is the position to 
be astrometrically measured by SIM, reflects the flux-weighted average of the two sources 
separately. Two classes of blending sources belong to this regime: a large fraction of binary 
companions to the source and all bright lenses. Of these two types of blending, Jeong, 
Han, & Park (1999) discussed in detail the effects of bright lenses on SIM observations and 
how to correct for blending effects caused by bright lenses. Here we focus on the effects of 
blending on events affected by binary star blending. 

When a source star is gravitationally amplified, it is split into two images located on 
the same and opposite sides of the lens, respectively (see Figure 2 of Paczyhski 1996). 
Due to the changes in position and amplification of the individual images caused by the 
lens-source transverse motion (see Figure 3 of Paczyhski 1996), the light centroid between 
the images changes its location during the event. The location of the image centroid relative 
to the source star is related to the lensing parameters by 

66 e =-^-(Tx + l3y), (4.2.1) 
tr + 2 

where T = (t — £o)Ae an d x and y are the unit vectors toward the directions which 
are parallel and normal to the lens-source transverse motion, respectively. If we let 
(x, y) = (S9 C>X , 59 C:V — b) and b = /3#e/2(/3 2 + 2) 1 / 2 , the coordinates are related by 

x 2 + V - 2 = a 2 , (4.2.2) 

where a = 9 E /2(/3 2 + 2) 1/2 and q = b/a = (3/((3 2 + 2) 1/2 . Therefore, during the event the 
trajectory of the source star image centroid traces out an astrometric ellipse (Walker 1995; 
Boden et al. 1998; Jeong et al. 1999). 

However, when an event is affected by binary star blending, the centroid shift is affected 
by the light from the companion star and its trajectory deviates from an ellipse. Therefore, 
by detecting the distortion of the astrometric shift trajectory, it is possible to identify the 
blend and thus correct for the blending effect. In Figure 3, we illustrate the distortion of 
the astrometric shift trajectory when an event is blended by a binary companion with a 
light fraction of / = 0.3. In the figure, the dotted line represents the unperturbed trajectory 
of the centroid shift with respect to the position of the lensed source star located at the 
origin. As mentioned, the trajectory is an ellipse. On the other hand, when the event is 
blended by a companion star, located at 56 B with respect to the lensed star, the position 
of the light centroid will shift toward the companion. In addition, the reference position of 
the astrometric measurements is not the position of the lensed source star but the center 
of light between the binary components before amplification. Due to these two effects, the 
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resulting trajectory (represented by a solid line) of the astro metric shifts is no longer an 
ellipse. Note that since the majority of binaries in this regime have long orbital periods, for 
example P ~ 30 yrs for a Galactic bulge binary system with A6 ~ 1 mas, the motion of the 
blend is not important. Assuming that the distribution of separations for binaries in the 
solar neighborhood determined by Duquennoy & Mayor (1991) can be applied to bulge and 
LMC populations, ^ 60% of Galactic and ^ 70% of LMC binaries belong to this regime. 

The trajectory of the centroid shift for a blended, binary event takes various forms 
depending on the fraction of blended light and the location of the blended star with respect 
to the lensed star. The centroid shift for an event blended by a binary companion is 
represented by 

^ c ,obs = ■ , I ~ f 77 Uabs^c - f59 B (A ahs " 1)1 • (4.2.3) 

J + ^absi-l - J) L J 

— * 

In the equation, the term including 59 c describes the elliptical motion of the centroid shift 
with respect to the lensed star (the elliptical term). On the other hand, the term including 
58b describes the linear shift caused by the light from the blended star (the linear term). 
Therefore, the shape of the astrometric shift trajectory for a blended event results from 
the combination of the elliptical displacement caused by gravitational lensing and the 
linear displacement toward the blended star. In Figure 4, we present various forms of the 
centroid shift trajectory as distorted by binary star blending. The left panels show how 
the trajectory changes from the unperturbed astrometric ellipse (in the top panel) with 
increasing fractions of blended light. To see the variation in the trajectory with respect 
to the location of the companion star, we also present the trajectories for various binary 
separations in the right panels. 



4.3. (Regime 3 X/d < A9 < 1.22X/D) 

In this case both the broad envelopes and fringes of the two PSFs are not matched 
(the lower right panels of Figure 1), and thus both centroids are shifted. The position of 
the centroid between the central broad envelopes is the flux-weighted average of the two 
sources and the amount of the shift is given by equation (3.2). One way to detect and 
correct for blending in this regime, therefore, is to measure the shift of the broad envelope 
centroid. However, the situation is slightly more complex for the position of the centroid of 
the narrow fringes. If the separation between two source is exactly the same as n times the 
width of the narrow fringe, i.e., A6 = nX/d, the 0th fringe of the blend coincide exactly 
with the nth fringe of the lensed source. However, this is true only at a specific wavelength. 
If the blended source is observed at twice as long a wavelength, for example, the 0th fringe 
of the blend coincides exactly with the 1st null of the lensed source. Therefore, from the 
difference in the interference patterns observed at multiple wavelengths, one can deduce the 
position and brightness of the blend even if the observer chooses to examine the central 
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fringe rather than the whole interference pattern. The sources of blending belonging to this 
regime include large-separation binary companions to the lensed source and very closely 
located field stars. 



5. Summary 

By observing gravitational microlensing events with SIM, one can correct for the effects 
of nearly all types of blending. With its diffraction-limited, <^ 0."5 imaging, a single mirror 
of SIM can resolve blended stars for a significant fraction of the events blended by field 
stars, such as regular and amplification-bias blended events. For events affected by very 
close blends with separations A9 ^ 10 mas, such as a large fraction of binary-star blended 
events, it is also possible to identify the blend by detecting a distortion in the astrometric 
shift trajectory. For events affected by blends with medium range separations, such as 
closely located field stars and wide binary companions, one can correct for the effects of 
blending by either detecting the shift of the broad envelope centroid between the PSFs of 
lensed and blended stars or by examining the narrow fringe patterns observed at multiple 
wavelengths. 

We would like to thank to M. Everett and P. Martini for a careful reading of the 
manuscript. We also would like to thank the referee (A. Gould) for useful suggestions that 
improved the paper. 
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Figure 1: The uncertainties in the lensing parameters for events whose high precision light curves 
are fit with models. The events are assumed to have £e = 15 days and (5 = 0.3, and are affected 
by blending with blended light factions of / = 0.3, 0.5, 0.7, and 0.9. The events are assumed to 
be observed 5 times/day during — 0.5£e < ^obs < 3£e with a photometric precision of p = 1%. The 
uncertainties in the lensing parameters are determined by computing x 2 , and the resulting x 2 as 
functions of tE/^E.o an d / are presented as contour maps. In each panel, the contours are drawn 
at la, 2a, and 3a levels from inside to outside. 
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Figure 2: The expected fringe patterns when a source star affected by a blended star with various 
separations is observed by SIM. To better show the centroid shift of the narrow fringe pattern, the 
region around the position of the central fringe of the lensed star is expanded in the lower panel 
of each PSF. For all of these model events, we assume a blended light fraction of / = 0.5 and 
arbitrarily normalize the intensity. 
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Figure 3: Illustration of the distortion of the astrometric shift trajectory by a binary companion. 
In the figure, the dotted line represents the trajectory of the centroid shifts with respect to the 
position of the source star located at the origin when the star is not affected by the binary-star 
blending. On the other hand, when the event is blended by a companion star, located at (0B,x,6B,y), 
the position of the light centroid will shift toward the companion. In addition, the reference position 
of the astrometric measurements is not the position of the lensed source star, but the center of light 
between the component stars. Due to these combined effects, the resulting trajectory (represented 
by a solid line) of the astrometric shifts is no longer an ellipse. This model event has lensing 
parameters of = 15 days, #e = 0.5 mas, and (5 = 0.5. The light fraction of the blended star is 
/ = 0.3. 
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Figure 4: Various forms of the centroid shift trajectory distorted by binary-star blending. The 
left-side panels show how the trajectory changes from the unperturbed astrometric ellipse (in the 
top panel) with increasing fractions of blended light. To see the variation of the trajectory with 
respect to the location of the companion star, we also show the trajectories for various binary 
separations in the right-hand panels. This model event has lensing parameters of t& = 15 days, 
#E = 0.5 mas, and f3 = 0.5. 



